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We present a framework for obtaining reliable solid-state charge and optical excitations and spectra from 
optimally-tuned range-separated hybrid density functional theory. The approach, which is fully couched within 
the formal framework of generalized Kohn-Sham theory, allows for accurate prediction of exciton binding 
energies. We demonstrate our approach through first principles calculations of one- and two-particle excitations 
in pentacene, a molecular semiconducting crystal, where our work is in excellent agreement with experiments 
and prior computations. We further show that with one adjustable parameter, set to produce an accurate bandgap, 
this method accurately predicts band structures and optical spectra of silicon and lithium flouride, prototypical 
covalent and ionic solids. Our findings indicate that for a broad range of extended bulk systems, this method 
may provide a computationally inexpensive alternative to many-body perturbation theory, opening the door to 
studies of materials of increasing size and complexity. 


Many solid-state systems exhibit strong excitonic effects, 
notably an optical excitation spectrum that is affected substan¬ 
tially by interaction between excited electron and hole quasi¬ 
particle states. The nature of this electron-hole, or excitonic, 
interaction is of central importance for a variety of applica¬ 
tions in, e.g., optoelectronics and photovoltaics ID. Neverthe¬ 
less, its accurate theoretical prediction remains a challenging 
task. It is common to account for such interactions within 
the framework of ab initio many-body perturbation theory, 
in which single-particle excitations are well-predicted from 
Dyson’s equation, typically solved within the GW approxi¬ 
mation Eia, and two-particle excitations are well-predicted 
using the Bethe-Salpeter equation (BSE) ||4]|5l. 

Current GW-BSE calculations are highly demanding and 
therefore presently impose significant practical limits on the 
calculated system size and complexity. Density functional 
theory (DET), in both its time-independent 013 and time- 
dependent (TDDET) lIMTD forms, is considerably more ef¬ 
ficient computationally. However, common (semi-)local ap¬ 
proximations to both DET and TDDET suffer from serious 
deficiencies which have precluded their use as a viable alter¬ 
native to GW-BSE in the prediction of excitonic properties 
II3. Eirst, quasi-hole and quasi-electron excitation energies 
are generally underestimated and overestimated, respectively, 
by the DET Kohn-Sham eigenvalue spectrum II3[T1. While 
the same functionals often perform better in the prediction of 
optical excitation energies of isolated molecular systems, the 
Kohn-Sham gap is typically similar to the optical gap HMH. 
In any case, they still fail in the solid-state limit ifTTl [12] [T9l - 
ED. Therefore, neither one- nor two-particle excitations are 
well-predicted in the solid-state, and hence the nature of exci- 
tons or their binding energies are not obtained. 

The failure of semi-local functionals in predicting solid- 


state absorption spectra has been traced back to an incorrect 
description of the long-range electron-electron and electron- 
hole interaction, manifested by the absence of a 1/q^ con¬ 
tribution 122] l23l to the interaction, where q is a wavevector 
in the periodic system. Several ingenious schemes for over¬ 
coming this deficiency have been suggested, including the use 
of an exchange-correlation kernel of the form — 

—Q;/(47r|r — r'|), where a is a system-dependent empirical 
parameter miEa; a static approximation to the exchange- 
correlation kernel based on a jellium-with-gap model ll26l : 
a “bootstrap” parameter-free kernel, achieved using self- 
consistent iterations of the random phase approximation 
(RPA) dielectric function EOl l27ll : a related ’’guided itera¬ 
tion” RPA-bootstrap kernel li28l : and the Nanoquanta ker¬ 
nel 1(121 l24l [29] l30l . derived by constructing the exchange- 
correlation kernel from an approximate solution to the BSE. 
Each correction provides a major step forward. However, 
none is a fully DET-based solution, as single quasiparticle ex¬ 
citations are obtained from GW, RPA, a DETh-U approach, or 
a scissors-shift correction. 

A different path for enabling TDDET calculations in the 
solid state is the use of (global or range-separated) hybrid 
functionals. These are still well within density functional 
theory, using the generalized Kohn-Sham (GKS) framework 
El (HIED, and their non-local Eock-like exchange com¬ 
ponent assists in the inclusion of long-range contributions. 
Although the time-dependent GKS equations have yet to be 
formally derived, hybrid functionals are already widely used 
for calculating optical properties. Eor gas-phase molecules, 
hybrid functionals can improve optical excitation energies, 
although standard hybrids still do not provide for accurate 
single-particle excitation energies EUTSlfTsl . TDDET using 
the Heyd-Scuseria-Ernzerhof (HSE) short-range hybrid func- 
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tional 1 ^ . where non-local exchange is introduced only in 
the short-range, can improve the absorption spectra of semi¬ 
conductors and insulators ll^ . although some discrepancies 
remain. However, the HSE functional still does not provide 
the desired long-range non-local contribution. The B3LYP 
hybrid functional ll34ll . in which a global 20% fraction of 
exact-exchange is used, was recently shown to yield TDDFT 
optical spectra for semiconductors in good agreement with ex¬ 
periment ll35l . but in some cases a larger fraction of exact ex¬ 
change was needed 061 . Although in this case a non-local 
contribution to the kernel tail does exist, it is global and pa¬ 
rameterized for a finite set of small organic molecules. How¬ 
ever, global and short-range hybrid functionals were shown 
to be insufficient predictors of band-structures in solid-state 
systems EtI, notably for molecular crystals OSl where ex- 
citonic effects are strong. Recently, Yang et al. 09l sug¬ 
gested a screened exact-exchange (SXX) approach, in which 
the local part of the hybrid calculation is set to zero and 
the time-dependent Hartree-Fock exchange is scaled down 
non-empirically per system by using the inverse of the di¬ 
electric constant, based on a ground state obtained from a 
scissor-corrected local density approximation (LDA) calcu¬ 
lation. Again, this led to improved performance for more 
strongly bound excitons. 

Ideally, we seek a DFT-based method where accurate one- 
and two-particle excitations can be read directly off of the 
eigenvalues of the time-independent (G)KS and the linear- 
response time-dependent (G)KS equations, respectively, us¬ 
ing a single exchange-correlation functional, from which a 
consistent exchange-correlation potential and kernel are de¬ 
rived. This challenge is not met by any of the above-surveyed 
methods. Recently, it was met for gas-phase systems, using 
the optimally-tuned range separated hybrid functional (OT- 
RSH) approach iflTl 14011411 . where the long- and short- range 
fraction of Fock-exchange is tuned non-empirically so as to 
obey rigorous physical constraints. This approach, elaborated 
further below, was shown to yield excellent fundamental and 
optical gaps for molecules ll42l l43l . More recently, it has 
been generalized so as to provide accurate single-particle ex¬ 
citations for both molecules Il44l - l46l and molecular crystals 
EH ill, and it was shown to capture gap renormalization in 
molecular solids ESl . Can this approach, then, resolve the 
long-standing challenge of providing an accurate one- and 
two-particle excitation spectrum in solid-state systems fully 
within the framework of (TD)DFT? 

In this Letter, we present a solid-state OT-RSH approach 
that achieves just that. It does so with an exchange-correlation 
potential and kernel that are fully consistent with the choice 
of the exchange-correlation energy, being its first and second 
functional derivative with respect to the density. We prove the 
accuracy of our approach by performing non-empirical cal¬ 
culations for pentacene, a prototypical molecular crystal and 
showing that it provides excellent agreement with GW-BSF 
calculations. Furthermore, with one empirical parameter - set 
to reproduce the known fundamental gap - we again achieve 
results that are comparable with both GW-BSF and experi- 
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FIG. 1. (Color online) (a) Band-structure (left) and density of states 
(right) of the pentacene solid, calculated using LDA (gray, dashed 
lines), GoWo@LDA (red, dashed lines), and OT-SRSH (black, solid 
lines). For all methods, the middle of the bandgap is shifted to 
zero, (b) The imaginary part of the dielectric function of the pan- 
tacene solid, with incident light polarization averaged over the a, b, 
and c main unit-cell axes, calculated using TDLDA (gray, dashed 
lines), Go Wo/BSE (red, dashed lines), and TD-OT-SRSH (black, 
solid lines). For visualization purposes, the leading absorption fea¬ 
ture (between 0.5 to 2.5 eV) was multiplied by a factor of 10 with all 
computational methods used. The OT-SRSH and TD-OT-SRSH re¬ 
sults were obtained using the parameters 7 = 0.16Bohr“^,a = 0.2, 
and e = 3.6. For computational details and convergence informa¬ 
tion, see the SI. 


ment for bulk silicon and LiF. The approach therefore emerges 
as promising for obtaining photoelectron and optical proper¬ 
ties accurately and efficiently for a broad range of extended 
systems. 

In the range-separated hybrid approach, the Coulomb inter¬ 
action is range-partitioned ll48l . e.g., via: m 

1 a-F/3erf(7r) ^ 1 — [a-F/?erf( 7 r)] 

— ‘ ' 4 “ ' * ' ( 1 ) 

r r r 

where a, /3 and 7 are parameters, and r is the inter-electron co¬ 
ordinate. The exchange expression corresponding to the first 
term on the right-hand side of Fq. 0 is then treated as in 
Hartree-Fock theory; the exchange expression corresponding 
to the second term is treated within the Kohn-Sham frame¬ 
work, typically using the LDA or the generalized gradient ap¬ 
proximation (GGA). In Eq. 7 is the range-separation pa- 
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rameter, i.e., it controls the range at which each of the terms 
dominates; a and f3 dictate the limiting behavior of the Fock- 
like exchange, which tends to a/r for r —>■ 0 and to (a + /3)/r 
for r —> oo. The resulting exchange-correlation energy is of 
the form; 

= (1 - a)E^l + aEl^ + (1 - (a + /3))i?^f,+ 

{a + I3)E^^ + Eksc, 

( 2 ) 


where KSx and KSc denote (semi-)local KS exchange and 
correlation respectively, and xx is a Fock-like exchange. 
SR and LR label short- and long-range terms, in which the 
Coulomb interaction is scaled using error functions. Within 
the GKS framework, the potential corresponding to the (semi- 
)local energy components is then obtained as a functional 
derivative, whereas the potential corresponding to the xx en¬ 
ergy components is obtained as a non-local Fock-like opera¬ 
tor. 

Here we generalize this approach to the time-dependent 
case, within the usual linear-response formalism of Casida ||9l 
l50l and the Tamm-Dancoff approximation. This is achieved 
by coupling GKS electron-hole pairs via an exchange- 
correlation kernel. This leads to an eigenvalue equation in 
which the eigenvalues are related to optical excitation ener¬ 
gies and the eigenvectors can be used to compute oscillator 
strengths, so that the complete optical absorption spectrum 
can be computed. The part of the TDDFT kernel originat¬ 
ing from the Hartree and (semi-)local Kohn-Sham potential in 
the ground-state DFT can be expressed as: 


a* 


Tj + (1 - -f (1 - a - 13) 


\r — r‘ 


LR 

xc 


\bj) (3) 


where = ^^S{r - r'), = §^(5(r - r'), and 

are the short- and long-range contributions of the 
(semi-)local KS exchange-correlation potential; and where 
a,b and i,j denote occupied and unoccupied states, respec¬ 
tively. The non-local exchange potential in the ground-state 
GKS leads to an additional term in the kernel, of the form 


- {ab\ 


erfc( 7 (|r — r'|) 


+ {a + 13) 


erf( 7 (|r-r'|) 


lu) 


(4) 

(see supplementary information (SI) for additional formal de¬ 
tails). 

As discussed above, optimal-tuning of the RSH parameters 
was shown to be crucial for achieving accurate description of 
molecular single-particle and optical excitations. This tun¬ 
ing procedure is challenging in the solid-state, as it involves 
a calculation of the system’s ionization potential and electron 
affinity from total energy differences, a problematic procedure 
for periodic systems (see m, and references therein). For 
the special case of molecular crystals, however, it was shown 
ll^ that predictive bandstructures can be achieved if 7 and 
a are optimally-tuned so as to obey the ionization potential 
theorem for an isolated molecule, with (3 chosen such that 


a + f3 = 1/eo, where Eq is the scalar dielectric constant, itself 
computed from first principles. We refer to this procedure as 
the optimally-tuned screened-KSH (OT-SRSH) approach. 

To test the efficacy of this approach for optical properties, 
we examine pentacene, a molecular semiconducting crystal 
of extreme interest in organic electronics and photovoltaics. 
For pentacene, the optimal tuning parameters are found to 
be a = 0.2, (3 = 0.08 (cotTesponding to Eq = 3.6), and 
7 = 0.16 Bohr~^ ll38]l . With these parameters we con¬ 
struct the appropriate time-independent (Eq. and time- 
dependent (Eqs. equations. These are then solved with 

the PARATEC ll5^ and BerkeleyGW ll5^ codes, which we 
modified to handle range-separated hybrids (see SI for a com¬ 
plete discussion). 

Eigure [T] (a) shows the resulting band structure and den¬ 
sity of states (DOS) calculated from EDA, OT-SRSH, and 
GW eigenvalues. The band gap predicted by OT-SRSH is 
1.9 eV, in good agreement with the GW gap of 2.1 eV. As 
expected, both improve drastically on the EDA gap, which is 
0.7 eV. Save for the slight difference in gap values, the OT- 
SRSH and GW band-structure and DOS are remarkably sim¬ 
ilar, demonstrating for this system that the eigenvalues of the 
OT-SRSH method are quantitatively useful approximations to 
single-quasiparticle excitation energies, as they are close to 
GW quasiparticle energies for several eV away from the band 
edges. Eigure [T] (b) presents the resulting imaginary part of 
the dielectric function, Im(E), calculated using TDLDA, TD- 
OT-SRSH and BSE, with the response averaged over incident 
light polarization along the a, b and c directions of the pen¬ 
tacene lattice (see SI for additional details). Importantly, our 
TD-OT-SRSH results do as well as GW-BSE in providing the 
expected picture of excitonic binding, i.e., comparison of the 
results with and without electron-hole interactions suggests an 
exciton binding energy of 0.45 eV from both methods. Eur- 
thermore, the BSE result agrees well with previously reported 
ones, e.g. those of Refs. EJED. The first singlet excita¬ 
tion is predicted to be at 1.46 eV in TD-OT-SRSH and 1.64 
eV in BSE. This small quantitative difference is essentially 
within the expected accuracy of either calculation. Moreover, 
this difference is consistent with that computed at the single¬ 
particle excitation level, and likely is inherited from it. In a 
similar manner, all presented optical excitation energies re¬ 
sulting from TD-OT-SRSH and BSE are very similar, with 
remaining differences within the desired accuracy. However, 
there are differences in the oscillator strength in part of the 
spectrum, notably at the first peak between 1.5 eV and 2.5 eV. 
These differences are primarily associated with the a-axis di¬ 
rection of incident light polarization. However, the scale of 
these differences is small when compared to the entire spec¬ 
trum. (Note that this spectral range is enhanced by a factor of 
10 in Eigjljb) to be visible.) Overall, then, the TD-OT-SRSH 
spectrum, while not identical to that of BSE, is in very good 
agreement with it. 

We now turn to two prototypical covalent and ionic bulk 
solids: Si and LiE, respectively. Eor these systems, we cannot 
use the original tuning procedure, because it was designed for 
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FIG. 2. (Color online) Top: The band-structure of bulk (a) silicon and (b) LiF, calculated using LDA (gray, dashed lines), OT-SRSH (black, 
solid lines), and GoWo@LDA (red, dashed lines). The Si bandstructure is also compared to empirical pseudopotential values taken from Ref. 
(21 (blue, dashed lines). For all methods, the middle of the bandgap is shifted to zero. Bottom: The imaginary part of e for (a) silicon and (b) 
LiF, calculated using TDLDA (gray, dashed line), TDSRSH (black, solid line) and GoWq/BSE (red, dashed line), and compared to experiment 
(blue, dashed line. Ref. 1551 for Si, and Ref. 1561 for LiF). The optimized SRSH and TDSRSH parameters are 7 = 0.11 Bohr“^, a = 0.2 , 
and e = 12 for silicon, and 7 = 0.58 Bohr'^, a = 0.2, and eq = 1.9 for LiF. For computational details, see SI. 


systems where inter-molecular hybridization is small. How¬ 
ever, it still possible to set a = 0.2 as a universally useful 
amount of short-range exact exchange Hanoi and demand 
a + /3 = I/eqj as before. We are then left with only one pa¬ 
rameter, the range-separation parameter, 7 . Here, we simply 
choose 7 so as to obtain the fundamental gap. We believe that 
future work may teach us how to obtain this parameter from 
first principles as well. 

As shown in the top panel of Fig|^a), the SRSH Si band- 
structure is fully comparable to the canonical empirical pseu¬ 
dopotential work of Chelikowsky and Cohen (54), as well as 
to GW results, whereas with LDA the gap, as expected, is 
too small. From the bottom panel we see that, in agreement 
with previous work Il24ll25]l . the TDLDA spectmm does not 
reproduce experiment satisfactorily. However, once again the 
TDSRSH result does almost as well as GW-BSE in predict¬ 
ing the experimental results, without any scissor shift or other 
correction operators at the single-quasiparticle level. Qualita¬ 
tively, the TDSRSH lineshape is indistinguishable from that 
of GW-BSE. Quantitatively, the peak positions are similar to 
within 0.05 eV. A similar picture emerges for LiE, as both 
single-particle and two-particle excitations are very close to 
those of GW/BSE. Specifically, both TDSRSH and BSE pre¬ 
dict the first and largest excitation at 12.6 eV, very close to the 
experimental value at 12.75 eV, and the overall spectral shape 
is satisfying. Here too, then, SRSH and TDSRSH are shown 
to capture excitonic effects, even though the Si excitons are 
much more delocalized and weakly bound compared to the 


case of pentacene, and the LiE excitons are strongly-bound 
and known to be highly challenging for many of the TDDET 
methods discussed in the introduction ll^ l26l l28l [39l . 

What are the physical origins of this success of the TD¬ 
SRSH approach in the solid state? Eirst, the ground-state 
calculation, being a generalized Kohn-Sham one, is capa¬ 
ble in principle of describing quasi-particle excitations owing 
to the non-local potential operator. Range-separation com¬ 
bined with long-range dielectric screening allows us to ful¬ 
fill asymptotic potential constraints while retaining the cru¬ 
cial balance of short-range exchange and correlation compo¬ 
nents, thereby making such prediction sufficiently accurate 
also in practice. Second, the asymptotic form of the TDSRSH 
kernel generates the desired non-local 1/q^ contribution by 
construction. Eurthermore, it is already scaled correctly via 
the non-empirical 1/eo parameter. Third, the non-locality of 
the TDSRSH kernel also alleviates the need for frequency- 
dependence, which would be necessary for bound exciton pre¬ 
diction with (semi-)local exchange-correlation kernels ll 20 l . 
Eor these reasons, highly accurate single-quasiparticle and 
two-quasiparticle properties of extended systems can be ob¬ 
tained in a predictive manner within the unifying framework 
of DET at computationally-modest cost. Last but not least, 
Onida et al. have already noted some time ago that “both 
the Green's functions and the TDDET approaches profit from 
mutual insight.” Here, we believe we achieve an important 
milestone towards that vision: on the one hand, the above 
work can be justified entirely from (TD)DET reasoning. On 
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the other hand, it is clear that the work has been motivated by 
the need to achieve the elegant quasi-particle picture obtained 
so naturally within many-body perturbation theory and that 
by achieving this goal we have created an effective simplified 
framework mimicking this picture. 

In conclusion, we have presented a new approach for quan¬ 
titative determination of single- and two-particle excitations 
in solids, based on range-separated hybrid density functional 
theory. The approach is fully couched within the formal 
framework of generalized Kohn-Sham theory. Furthermore, 
it is based on a fully self-consistent choice of exchange- 
correlation energy, potential, and kernel. We have shown that, 
with one empirical parameter at most, it produces results that 
are in quantitative agreement with those of many-body pertur¬ 
bation theory in the GW-BSE approximation, for three proto¬ 
typical systems: Si, a covalent solid, LiF, an ionic solid, and 
pentacene, a molecular solid. In particular, it fully captures 
excitonic interactions for both strongly and weakly bound ex- 
citons. In doing that, it answers a long-standing challenge of 
(TD)DFT - providing an accurate one- and two-particle ex¬ 
citation spectrum for solid-state systems within a fully con¬ 
sistent framework. We envision that it could emerge as a 
useful low-cost substitute to many-body perturbation theory, 
as well as provoke further developments within density func¬ 
tional theory. 
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